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We present a quantum theoretical treatment of light-matter coupling in the system consisting of 
a quantum dot and a spherical core-shell metal-dielectric multilayer nanoparticle. It is shown that 
both weak and strong coupling regimes can be realized in the set-up. Specifically, we demonstrate a 
strong coupling regime between a quantum dot and a nanoparticle, when the quantum dot resonance 
is tuned to the frequency at which normal component of effective nanoparticle permittivity is crossing 
zero. Moreover, we demonstrate the regime at which the quantum dot decays much faster than in 
vacuum (due to the large Purcell factor) and at the same time radiates more power to the far field. 

This findings pave the way towards more efficient control over radiation properties of quantum 
emitters. 


Tailoring the radiation properties of the quantum emit¬ 
ters by placing them inside or in the vicinity of the micro¬ 
cavities has recently emerged in a rapidly evolving field 
known as Cavity Quantum Electrodynamics^ - — . Typi¬ 
cal solid state realizations of microcavities are based on 
photonic crystals with defect (s),— or Bragg mirrors^. To 
boost the light-matter interaction in the sub-wavelength 
regime, there has been some on-going interest in the re¬ 
alization of strong light-matter coupling with plasmonic 
cavities that are well-known for confining light below the 
diffraction limit (see Ref— i, for a comprehensive review). 

Depending on the degree of coupling between the emit¬ 
ter and the external medium, the situation can be clas¬ 
sified into ‘weak-coupling’ and ‘strong-coupling’. In the 
former case, energy propagation from a quantum emitter 
to it’s environment is unidirectional, i.e. energy of the 
emitter decays in a monotonous way, while in the latter 
case, the emitted energy is reflected back to the original 
emitter and then re-emitted by the emitter and so on, 
thus the energy transfer is bi-directional. The QE can 
be a quantum dot (QD), fluorescent molecule, or a sin¬ 
gle atom. In the study of weak-coupling, the external 
body near the emitter is considered as perturbation to 
the original system - atom-in-vacuum. However, to un¬ 
derstand the full physics in the strong-coupling regime, 
non-perturbative calculations are imperative. To date, 
strong-and ultra-strong-coupling have been studied for 
wide variety of systems^ - ™. In this Communication, we 
consider the coupling of a quantum emitter to a multi¬ 
layered spherical core-shell nanoparticles. 

Core-shell nanoparticles^S have recently become an 
extremely popular object of study in nanophotonics. 
Namely, core-shell nanoparticles have been suggested for 
the realization of both super-scattering and total scat¬ 
tering cancellation^. Moreover, it has been shown that 
the multilayered metal-dielectric core-shell nanoparticles, 
can be used for the realization of spherical hyperbolic 
cavities, i.e. cavities made from hyperbolic metamateri- 
als22. Hyperbolic metamaterials are a particular class of 
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metamaterials where the dielectric tensor has indefinite 
signature, and the iso-frequency surface has infinite vol¬ 
ume that leads to broadband singularity of photonic local 
density of states (LDOS), which results in the broadband 
increase of the Purcell factor for the emitters planced 
in the vicinity of hyperbolic metamaterials^ - —. At the 
same time for the case of hyperbolic metamaterials, the 
enhanced Purcell factor does not lead to the enhance¬ 
ment of the radiation intensity of the emitters placed 
in its vicinity. Contrary to that, radiation intensity is 
actually decreased, since most of the emitter energy is 
transferred to the near-field modes of the metamaterial, 
which decay non-radiatively. In this Communication we 
show that substitution of the bulk hyperbolic metamate¬ 
rial with a finite-size hyperbolic cavity can provide both 
large Purcell factor and radiation intensity increase for 
the quantum emitters placed in its vicinity. 


Hyperbolic Metamaterial 
Cavity 



FIG. 1. A hyperbolic cavity of spherical shape (radii a) is 
in the vicinity of a quantum dot at a distance h from the 
surface of the sphere. The far-held spectrum is attained via 
a detector at distance rn from the center of the sphere. The 
above schematic shows a ‘cut-view’ of the sphere. 

While a considerable number of papers have been ded¬ 
icated to the studies of the Purcell effect in hyperbolic 
metamaterials, the model of the strong coupling of the 
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quantum emitter to the hyperbolic cavity has not been 
put forward so far as to our knowledge. 

In this Communication, we report that a hyperbolic 
cavity and a quantum emitter exchange photon in the 
strong-coupling regime. Via non-perturbative calcula¬ 
tion, we reveal the existence of two distinct anti-crossings 
in the local polarization spectrum: one at the vicinity of 
the nanoparticle pseudomode frequency formed by the 
series of high angular momentum modes, and one in the 
vicinity of the frequency, at which the metal-dielectric 
cavity experiences the topological transition - from the 
elliptical to hyperbolic regime. However only in the latter 
case, the radiation reaches the far field zone. Moreover, 
the mode that is detectable at far-held has quite large 
Purcell factor. We report an unique behavior of the 
strong-coupled system - radiation power enhancement 
and large Purcell factor at the same spectral position in 
hyperbolic metamaterial cavity. 

In order to understand the radiation properties of the 
quantum emitter to the environment around it as shown 
in FigCD the Purcell factor T/ro and Lamb shift Aw are 
calculated upto the leading order of perturbation which 
are shown in Fig. [2l for two different orientations of quan¬ 
tum dot. The quantum emitter is a point quantum dot 
and is placed at a distance h from the surface of the spher¬ 
ical cavity made of alternating dielectric and metal con¬ 
centric layers of fixed filling fraction and the following pa¬ 
rameters are fixed throughout this paper unless otherwise 
mentioned: the sphere consists of ten periods and the 
permittivity of the dielectric is £\ = 2.25 and the metal 

UJ 2 

is modeled by Lorentz model, e(w) = 1-1- 2 —; 

the parameters for the model are - w p /wt = 0.5 and 
y/w t = 10 -4 and h = 10 -3 x A t , where A t = 27rc/w t . The 
radius of the multilayered particle is 0.1 x At, the filling 
fraction of metal in one period is 0.5, and the linewidth 
of the quantum dot is 10 -4 x w*. 

Purcell factor can be calculated via 


r = r n 


2w 21_c- 


HeqC 2 


d • G(r QD ,r Q D,W2i) - d 


(1) 


where To is spontaneous decay rate in vacuum, W 21 is 
quantum dot frequency, W21 plus the Lamb shift, Aw. In 
the leading order of perturbation, W 21 can be replaced by 
‘bare QD frequency’, W 21 and d is dipole matrix element 
of QD while G(r,rQ D ;w) is the Green’s function of the 
system - 


V x V x — 



G(r,r QD ;w) =T( 5 (r-r QD ), ( 2 ) 


Since the quantum dot has been placed at a very close 
proximity to the spherical cavity (compared to the wave¬ 
lengths of interest), a number of modes is excited by 
the QD that give rise to a series of narrow peaks in ad¬ 
dition to the primary peaks at W 2 i/wt « 1.04andl.l2. 
The second dominant peak (w 2 i/wt ~ 1.12) in the Pur¬ 
cell factor and Lamb shift can be understood via the ef¬ 
fective medium theory for metamaterial homogenization 
that reads - 


— d\ Em T ^2 Ed : / a\ 

e- 1 = di£^ +d 2 £j 1 - 

where e r is the radial component of the dielectric ten¬ 
sor and Eg :V the transverse component of the dielectric 
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I being the usual 3x3 unit dyad. The Green’s function 
in Eq. [2] is calculated by expansion of the field as a series 
of spherical harmonics^. 

The Lamb shift reads as 


SR 


Aw = — - 


d • G(r QD , r QD ; w) • d 


UJ ; 


21 


fen 


(3) 


FIG. 2. Purcell factor and Lamb shift for a dipole along (a) x- 
directed dipole, (b) ^-directed dipole, for a hyperbolic cavity 
shown in Fig. [T] Since, the quantum dot has been placed at 
a very close proximity to the sphere, a number of modes is 
excited by the QD that give rise to a series of narrow peaks 
along with the peak at epsilon-near-zero frequency (W 21 /wt = 
1.118) and at W 2 i/wt = 1.038. The shaded green regions 
represent the hyperbolic regime. Note the logarithmic scale 
of the left axis for Purcell factor. 
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tensor. Two regions can be distinguished from the above 
description, namely 1) elliptical regime where e r and £g tV> 
have same sign, 2) hyperbolic regime where e r and eg tV 
have opposite signs. The hyperbolic regime is shaded 
by green in Fig. [2j the unshaded portion represents the 
elliptical regime. The enhanced Purcell factor around 
u>/uj)t « 1.12 occur at the the spectral point where hy¬ 
perbolic to elliptical topological transition occurs. As will 
be discussed later that the first resonant peak originates 
from the hybridization of higher order spherical harmon¬ 
ics. The Lamb shift for the x and ^-oriented dipole mo¬ 
mentum is rather large compared to typical Lamb shift 
in experiment.— This high Purcell factor and large Lamb 
shift is the initial signature of strong-coupling that re¬ 
quires further scrutiny via non-perturbative calculation. 

Before embarking on the non-perturbative calculation, 
we would like to mention one little subtlety with the 
real part of the Green’s function due to it’s diverging 
nature. However, this divergence can be gotten rid of 
by setting high-/c cut-of£I. This regularization extracts 
realistic Lamb shift because any finite-sized dipole must 
have a finite Lamb shifts. Moreover, near the resonance, 
the homogeneous vacuum contribution is orders of mag¬ 
nitude smaller than the scattering part of the Green’s 
function. Thus disregarding the homogeneous part of 
Green’s function is a reasonable approximation and this 
allows one to get rid of the unphysical divergence of Lamb 
shift. In fact, if the whole calculation is performed us¬ 
ing Lippmann-Schwinger approach, the divergent part of 
Green’s function is actually gotten rid of and it is the 
scattering part that contributes eventually^. Hence, in¬ 
stead of total Green’s function G(r, r; cj ) that consists of 

vacuum Green’s finction G (r, r; u>) and the scattering 

— sc 

part of the Green’s function G (r, r; w), it is reasonable 
to replace G(r, r; oj) by G (r, r; w) in Eq. [T| and [3] 

Now, to study the coupling properties in a non- 
perturbative fashion, let us consider a state-space T-L with 

two states^2 namely, the ground state, |1) = ^ ^ } ’ and 

the excited state, |2) = 

Considering weak excitation, i.e. only one quantum 
of light is associated at a time, the Hamiltonian of the 
system can be written as — 



E / dr / dw; fuxi f A (r,wj) fA(r,Wi) 

A=e,m 

- [cr + d + cr - d] • F(r QD ) + foj2icr + <? - , 


where e and m stand for electric and magnetic field, re¬ 
spectively, that are related to the continuum Bosonic- 
field operator f(r,w;) and it’s Hermitian conjugate 

f\r,u;;) and, with eigenfrequency oji. The operators 
d + = |2}(1| and a~ = |1)(2| are the usual Pauli oper¬ 
ators for exciton, and d being the dipole matrix element 
(i.e., d = (2|er|l)) between the ground state, |1) and 
excited state, |2) with transition frequency, ui 2 i- Since, 
the primary interest of this work is in the optical regime, 


we set the relative permeability, /r = 1 and exploiting 
bosonic operators, the field operator takes the following 
form - 

F(r) = E°(r) + — • ^-G(r, r QD ; w) • d 

E 0 c 2 (5) 

x [ct + (w) +a~(u})] , 

where E (r) is free-held operator, i.e. the held without 
the presence of quantum emitter. The local polarization 
spectrum that provides insight on quenching and other 
propagation mechanisms at the dipole position tqd can 
be expressed via - 


P(r QD ,w) = (ct + (w)ct (u)) 


W QD - w 2 - ^d • G(r, r QD ; w) • d/fe 0 


( 6 ) 


In Fig. m, a) and (b), the local polarization spectrum 
of the quantum dot are shown as a function of transition 
frequency of quantum dot, wqd- The ‘red-stars’ repre¬ 
sent spectral position of the quantum dot. As seen from 
Fig. [3] (a) and (b), two distinct anticrossings are evident, 
at two different spectral regimes. Note that the anti¬ 
crossings in Fig. [3Ka) and (b) are centered around the 
two spectral points with high Purcell factor. This is the 
direct signature of strong-coupling of photon and HM. 
Despite the strong-coupling (i.e. strong light-matter en¬ 
ergy exchange), the resonance at topological transition 
(lo/oj t 1.12) is rather broad. This fact seems counter¬ 
intuitive. The resonance at the second anti-crossing has 
Fano line-shape: a slow varying broadband resonance in¬ 
teract with a narrow resonance that results the broad 
asymmetric line-shape. This explains the broad nature 
of the line-shape despite strong-coupling. 

To understand the distinctive nature of the two anti¬ 
crossings the radiative property is further explored. The 
non-local spectrum of the QD, i.e., spectral density at a 
point r(^ tqd) for a QD position at tqd reads^, 


S( r,w) 



d t 2 (F f (r, £i)F(r, t 2 )) 


, e *“(t 2 —ti) 


( 7 ) 


Exploiting the field operator F(r, t\) of Eq.0 the far-held 
spectrum in Eq. [7] reads, 

c , n ' G(r, tqd; w)(w + wqd)/eo 

S(r, yqd,ui) = - 2 -5--- 

W QD - w 2 - ^j-d • G(r, r QD ; w) ■ d/fe 0 

( 8 ) 

Eq. [8] offers an intuitive picture since it incorporates 
the propagator G(r, tqd;w), that delineates the photon 
propagation to any point r (^ tqd) from the quantum 
emitter at rQD- 

The far-held spectrum at Fig. m c) and (d) provides 
interesting insight all the modes that appear in the 
near-held polarization spectrum in Fig. [3] (a) and (b), do 
not propagate to the far-held. The mode near topologi¬ 
cal transition (i.e., W2i/wt~1.12) is the only one that has 
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FIG. 3. (a) and (b), Local polarization spectrum of a quantum dot for different spectral position of quantum dot. (c) and 
(d) The far-held spectrum for the same spectral position of quantum dot same as (a) and (b). As shown here, the modes at 
ui/u>t « 1.04 does not appear at far-held while the mode near the topological transition of hyperbolic metamaterial. 


far-field signature. Since, the modes at first anticross¬ 
ing (co/uJt ~ 1.04) are whispering gallery modest that 
do not propagate to the far-field, and they show zero 
spectral signature of this mode at far-held. However, the 
mode at topological transition provides large Purcell fac¬ 
tor while it transfers photon to the far-held. Unlike the 
traditional HMMs, this demonstrates large Purcell fac¬ 
tor accompanied by radiation power enhancement. This 
is the central result of this work. 

To further substantiate the aforementioned claim, har¬ 
monic decomposition of the Green’s function in spherical 
harmonics has been performed (results not shown here). 
It is found that the hrst anticrossing has no dipole con¬ 
tribution (i.e., the contribution of J = 1 is negligible), 


contrary to the mode at topological transition that has 
significant dipole contribution in the Green’s function. 
We refer the mode at W 2 i/u>t ~ 1.04, is a pseudomode 
that appears as a contribution from higher order har¬ 
monics, similar to Ref4. 

To summarize, strong light-matter interaction in a hy¬ 
perbolic cavity is demonstrated where two distinct an¬ 
ticrossings are present in the local polarization spec¬ 
trum. The mechanisms responsible for the appearance 
of the two anticrossings are distinctive - the first anti¬ 
crossing (i.e., the pseudo-mode) appears as contribution 
from higher order spherical harmonics of the cavity that 
is similar to the behaviour of plasmonic nano-particle, 
whilst the second anti-crossing appears due to the hy- 
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perbolicity of the hyperbolic cavity, at the topological 
transition of the cavity. The mode at topological transi¬ 
tion propagates photon to the far-held as opposed to the 
pseudo-mode. This is quite unlike the case of quantum 
dots in the vicinity of planar layered hyperbolic meta¬ 
materials, where the large Purcell factors co-exist with 
the lowered radiation of the far held. The extension 
of the analysis for multiple quantum emitters is quite 
straightforward via Dyson’s equatio n 4 ' 35 ! 36 that involves 


the dressing of Green’s function via incorporating the 
scattering by different quantum emitters and the rest of 
the system. 
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